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We report measurements of transfer functions and flux shifts of 20 on-chip high Tc DC SQUIDs 
half of which were made purposely geometrically asymmetric. All of these SQUIDs were fabri- 
cated using standard high Tc thin film technology and they were single layer ones, having 140 nm 
thickness of YBa2Cu307_a; film deposited by laser ablation onto MgO bicrystal substrates with 24° 
misorientation angle. For every SQUID the parameters of its intrinsic asymmetry, i. e., the density 
of critical current and resistivity of every junction, were measured directly and independently. We 
showed that the main reason for the on-chip spreading of SQUIDs' voltage-current and voltage-flux 
characteristics was the intrinsic asymmetry. We found that for SQUIDs with a relative large induc- 
tance (L > 120 pH) both the voltage modulation and the transfer function were not very sensitive 
to the junctions asymmetry, whereas SQUIDs with smaller inductance (L ~ 65 — 75 pH) were more 
sensitive. The results obtained in the paper are important for the implementation in the sensitive 
instruments based on high Tc SQUID arrays and gratings. 

PACS numbers: 74.50.+r 



I. INTRODUCTION 

The high-transition-temperaturc superconducting 
quantum interference devices (high-Tc DC SQUIDs) 
each consisting of two bicrystal Joscphson junctions are 
key elements for many sensitive instruments such as 
extremely low noise magnetometers^^, the series- and 
parallel- SQUID arrays^, the superconducting quantum 
interference grating^, etc. 

However, the further developments of these devices are 
limited by significant on chip spreading in the critical 
current and normal resistance of high Tc Josephson junc- 
tions which seems to be unavoidable for grain-boundary 
junctions4&±^H. 

The different values of critical current and normal re- 
sistance of two Josephson junctions result in turn in the 
vast spreading of the output voltage-current (VCC) and 
voltage flux (VFC) characteristics for on chip high Tc 
DC SQUIDsHiii^. 

The influence of the junction asymmetry on the output 
characteristics of high Tc DC SQUIDs has been analyzed 
i n 16 i 17 ' 18 i 19 ' 20 , where it was shown that the transfer func- 
tion for asymmetric SQUIDs can be substantially differ 
from that of symmetric SQUIDs. 

The extensive comparison of experimental characteris- 
tics of intentionally fabricated asymmetric high Tc DC 
SQUID with computer simulations has been performed 
ini 7 ., where in most cases a reduction of the voltage-to- 
flux transfer function of the asymmetric SQUID as com- 
pared to the symmetric SQUID has been observed. 

From the other point the output characteristics of high 
Tc DC SQUIDs can be used to infer the information 
about the on chip distribution of the critical currents 
and normal resistances of the Josephson junction o 11 ' 12 . 

In known experiment er 11 ! 12 ' 17 the asymmetry of critical 
current of two Josephson junctions in DC SQUID loop 
has been determined indirectly from the measured flux 



shift of the voltage-flux curve. Since for high Tc Joseph- 
son junctions the density of critical current and the resis- 
tivity of the junction are interrelated, it also influences 
the value of the last quantity. In addition, as is shown in 
the paper, this method has a restricted range of validity. 
That is why, as was noted ini 7 -, a clear-cut comparison 
between simulation and experimental data requires the 
independent experimental determination of the current 
and resistance asymmetry. 

Therefore, the main purpose of our paper is the in- 
dependent and direct determination of the current and 
resistance asymmetry and the investigation of the influ- 
ence of the junctions asymmetry on the scattering over 
the chip of the main output parameters of DC SQUID, 
i.e., its voltage modulation, voltage-to-flux transfer func- 
tion, flux shift. 

To this end we investigated 20 purposely made asym- 
metric high Tc DC SQUIDs. For the first time we mea- 
sured independently and directly the current and resis- 
tance asymmetry for every junction in the interferometer 
loop. We show that the main reason for the junctions 
asymmetry is the on-chip spreading in the resistivity and 
the critical current density. Even for a geometrically 
symmetric design of the interferometer the critical cur- 
rent densities of the two junctions may be substantially 
different. We found that, in general, both the voltage 
modulation and the transfer function are not very sensi- 
tive to the junctions asymmetry, whereas the flux shift 
of the VFC shows approximately linear dependance on 
asymmetry parameters. 

The paper is organized as follows. In Section II the 
current and resistance asymmetry for the junctions of 
DC SQUID are defined and the analytical expressions 
for transfer function and flux shift for asymmetric DC 
SQUID are given. The experimental part of the paper is 
in detail described in Section III. Section IV is devoted 
to an extensive analysis of the influence of SQUID asym- 
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FIG. 1: The asymmetric DC SQUID 



metry on the transfer function and flux shift of VFC. 
Obtained results are shortly summarized in Section V. 



II. ASYMMETRIC DC SQUID 

A. The asymmetry parameters 

An asymmetric DC SQUID shown in Fig[T] consists 
of a superconducting loop of inductance L intersected 
by two Josephson junctions, which have different critical 
currents lex, Ic2 and normal resistances i?2- The 
capacitances of the two junctions, which are not shown 
in FigQ] are assumed to be equal. In order to describe 
the junctions asymmetry we define the average values of 
the critical current Ic and resistance R, and a current 
asymmetry 7 and a resistance asymmetry p as follows: 

Ici = (1 + 7)Ic, I C 2 = (1 - i)Ic, Ri = + p), 

i?2 = R/(l — p), where 
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Note that the SQUID critical current Isq and its resis- 
tance Rsq are being measured directly from the voltage- 
current characteristic of the DC SQUID. 

There are two different origins of the junction asymme- 
try: geometrical asymmetry and intrinsic asymmetryi 7 -. 
For bicrystal grain boundary junctions geometric asym- 
metry is associated with the different width w of the junc- 
tions. We describe geometric asymmetry by the parame- 
ter a g according to W\ = (1 + a g )w and w% = (1 — a g )w, 
with w — (wi + u>2)/2. The intrinsic asymmetry is asso- 
ciated with different values of the current density jo and 
the resistivity po for two junctions. We describe intrinsic 
asymmetry by the parameters a.j and a p according tc 17 : 



3oi=3o0- + otj), Pi = Po/(l + OL p ) 

3o2=3oO--0!j), P2 = Po/(l - Dtp) (3) 

The parameters of intrinsic asymmetry ay and a p can 
be obtained from the independent measurements of the 
parameters of the bulk asymmetry 7, p and geometric 
asymmetry a^i: 
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Since high Tc Josephson junctions obey the scaling law 

IcR ~ jc^ 2 —, the parameters aj and a p are interrelated 
to each other— 
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B. Output voltage across asymmetric SQUID 

In general, the output voltage across asymmetric 
SQUID is the complicated function of several dimen- 
sionless parameters^: V = F(i,a,P,T,ipx, p,j), where 
i = I / Ic (I is the bias current); a = L/Lf {Lf — 
(<$>o/2Tty /ksT is the fluctuation inductance which is 
equal approximately to fOO pH at T=77 K); (3 = 
2LIc/$o; r = ZnksT / QqIc is the noise parameter 
(ks is Boltzmann constant, T is absolute temperature); 
ipx — tt&x/Qq {&x is the external magnetic flux, $0 
is the flux quantum). Three parameters a,/3,r are not 
independent, but are subject to the relation a = n(3T. 

There are several properties of the output characteris- 
tics of asymmetric SQUIDs which allow one to identify 
them by experiment 2 —. First, in the presence of applied 
flux the voltage across asymmetric SQUIDs is not an odd 
function of the bias current, i. e., V(—i) ^ — V(i). Thus, 
the quantity V(—i) + V(i) can be used as one of the mea- 
sure of SQUID asymmetry. It is a periodic function of 
the applied flux and it depends on the bias current and 
on the parameters of asymmetry. 

The second property which is known for a long time 
(see, for example, Refl221) is the shift of the voltage-flux 
characteristic (VFC) under reversal of the bias current. 
The shift is commonly attributed to the asymmetry in 
the critical current of the two junctions: 



A$ = L(Ici - I C 2) = 7/3*o 



(0) 



This property has been widely used for experimental de- 
termination of the current asymmetr y 11 ' 12 ! 17 ' 23 . How- 
ever, the expression ^ cannot describe some experimen- 
tal facts. In particular, it does not depend on the bias 
current while the experimental flux shift does depend on 
i. For large inductance DC SQUIDs the picture is more 
complicated^. There is no simple relation between cur- 
rent asymmetry and a flux shift. In addition, asymmetry 
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in resistance contributes also to the total flux shift. The 
expression for the flux shift under a reversal of the bias 
current is as follows 2 ^: 

l Q ,. r , , pa . 

-^_ = _e( J ,r,p, 7 ) + — , (7) 

The explicit expression for the quantity which is too 
cumbersome one can find in Ref. [20| (Eq. 31b). 

As is seen from ([7]) the flux shift is a complicated func- 
tion of the bias current and the parameters of asymme- 
try. The dependance on SQUID inductance L is hidden 
in the second term in ([7|). Due to the relation a — irf3T 
this term can be written as pi/3/2, which is equivalent to 
the expression found for the flux shift in the limit of a 
large bias current 

It is worth noting that for i > the voltage-flux curve 
is shifted to the negative side of the flux (A&x < 0) 
relative to a symmetric voltage-flux curve, while for the 
reversed bias (i < 0) the voltage-flux curve is shifted to 
the positive side of the flux (A<E>x > 0). 

C. Transfer function of asymmetric SQUID 

As is known, the shape of experimental VFC is not 
symmetric. The asymmetry of the junctions induces dis- 
tortion of the V($)- curves, which leads to different values 
of V£ and for the maximum positive and negative 
slope of the V(<&)-curves, respectively. The impact of the 
asymmetry on the transfer function has been studied by 
computer simulations in Ref. [T3- It was shown that the 
normalized transfer function — V$&q/RIc can be 
factorized as v§ = g ± {Yf3)f ± {^f 1 l3). The explicit analyt- 
ical expressions for the quantities g + iTfi) and / + (7,/3) 
for limited range of parameters T and (3 (T(3 < 1, (3 < 5) 
have been given in Ref. [I?]. 

III. EXPERIMENTAL 

1. High-Tc DC SQUID preparation and characterization 

The principal layout of an asymmetric high-Tc DC 
SQUID is shown in Fig. O It consists of the interfer- 
ometer loop and two Josephson junctions with different 
widths. A total of 20 such SQUIDs was integrated on one 
chip. All have the same superconducting loop width of 
8 /xm and slit width of 4 /im, respectively. Three differ- 
ent slit lengths of 25 /Ltm, 75 fim and 125 /jm were used, 
what partitions the SQUIDs into three groups: small (S 
SQUID), medium (M SQUID) and large ones (L SQUID). 
The slit lengths are reflected in SQUID inductance ranges 
of L-65-75 pH, L-125-135 pH and L-185-195 pH for the 
S, M and L SQUIDs, respectively. Within the groups, pa- 
rameters of asymmetry were varied, using various junc- 
tion widths w\ and 102 between 0.4 and 1.2 /xm. All 
SQUIDs were connected via a common superconducting 
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FIG. 2: Design of asymmetric high-Tc DC SQUID 




FIG. 3: Color online. Photography of three asymmetric high- 
Tc DC SQUIDs. 



line, which was used to inject a DC current Ii n j into the 
SQUIDs. Figure [3] shows a photograph of three SQUIDs 
connected that way. 

The SQUIDs were fabricated using a standard thin- 
film technology, which was described in detail in Ref. [2f| 
Using laser ablation, a 140 nm thick YBCO film is de- 
posited onto a 24° MgO bicrystal substrate with 10 mm x 
10 mm x 1 mm dimension. On top of the non-structured 
YBCO film a 100 nm thick gold layer was deposited by 
thermal evaporation and subsequently structured by lift- 
off. Finally, the YBCO layer was patterned using stan- 
dard e-beam lithography and ion-beam-etching. As al- 
ready shown in^, Josephson junctions with widths down 
to 0.4 /zm show no degradation in critical current den- 
sity. In addition, the chip was covered by a Teflon layer 
which protected it from water during thermo cycling 26 . 
Then the SQUID chip was glued on a non-metallic holder 
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FIG. 4: Color online. Determination of critical current from 
experimental VCC (thick line). The critical current Ic to 
which the experimental VCC is normalized, is varied until 
the normalized experimental VCC fits to one of the theoretical 
VCCs (thin lines). 



with a wire-wound coil beneath, which was used to apply 
magnetic flux to the SQUIDs. 

All measurements were performed in liquid nitrogen 
at 77 K in a well shielded environment. First, for ev- 
ery SQUID we have measured its voltage-current (VCC) 
and voltage-flux (VFC) characteristics. Second, one of 
the two junctions in the loop was measured directly. For 
that after all the measurements at the SQUIDs have been 
performed, all SQUID loops have been cut to allow the 
independent measurements of the characteristics of one 
junction per SQUID. To do that the right superconduct- 
ing arm of every SQUID loop was removed by additional 
optical lithography and chemical etching. This allowed to 
measure the critical current and normal resistance of the 
Josephson junction in the remaining arm of the SQUID. 
Therefore, we determined for the first time independently 
and directly the current density and resistivity for every 
junction in the interferometer loop which allowed us to 
investigate their influence on the SQUID characteristics. 

2. SQUID characteristics 

The first group of measurements concerned the char- 
acterization and classification of the SQUIDs. All these 
parameters are listed in Table [U 

Within the groups (S, M, L) the junction sizes wi and 
W2 are varied between 0.4 and 1.2 /im, yielding the ge- 
ometrical asymmetry a g . Here, the junction widths are 
layout parameters. The resulting SQUID inductance L 
was actually determined, however. This was done in the 
following way. 

The injection current Ij n j flows around the upper part 
of the SQUID loop (see Fig. [2]). This coupling part of the 
SQUID has an inductance of L c which can be measured, 
because the injection current produces a magnetic flux 
<!> C =L C li n j in the SQUID loop. The measured voltage- 



TABLE I: DC SQUID parameters; wi, W2 are junctions' 
widths; a g accounts for geometric asymmetry; L, Isq, Rsq 
are the DC SQUID inductance, its critical current and normal 
resistance; AVmax is the maximum swing of the VFC; l m ax 
is the bias current which provides AVmax- 



SQUID#group 


Wi,W2 

/j,m 


a g 


L 

pH 


Isq 
LiA 


Rsq 
Q 




I max 

/xA 


1 S 


0.6, 0.6 





71.4 


36 


3.95 


35.8 


29.5 


2 S 


0.6, 0.6 





71.4 


22 


5.79 


48.2 


22.1 


3 S 


0.8, 0.4 


0.33 


73.8 


35 


2.99 


21.8 


35.3 


4 S 


0.8, 0.8 





65.4 


29.5 


4 


38.8 


25.6 


5 S 


1, 0.6 


0.25 


66.6 


25.5 


4.82 


43.8 


23.1 


6 S 


1.2, 0.4 


0.5 


70.6 


42 


3.4 


35.6 


36.1 


7 M 


0.6, 0.6 





131.6 


31 


4.93 


21.5 


28 


8 M 


0.6, 0.6 





131.6 


33 


3.87 


16.4 


26.3 


9 M 


0.8, 0.4 


0.33 


134.1 


40 


4.56 


16.4 


32.1 


10 M 


0.8, 0.8 





125.7 


30.5 


4.11 


16.9 


25 


II M 


0.8, 0.8 





125.7 


56.5 


2.61 


12.2 


51.2 


12 M 


1, 0.6 


0.25 


126.8 


55 


3.1 


15.1 


45.1 


13 M 


1, 0.6 


0.25 


126.8 


35 


3.27 


14.5 


29.7 


14 M 


1.2, 0.4 


0.5 


130.5 


42 


3.55 


15.5 


39.5 


15 M 


1.2, 0.4 


0.5 


130.5 


53 


2.76 


13.7 


48.3 


16 L 


0.6, 0.6 





191.5 


29 


4.34 


9.2 


24.2 


17 L 


0.8, 0.4 


0.33 


194.1 


25.5 


4.68 


8 


19 


18 L 


0.8, 0.8 





185.8 


46 


3.41 


6.7 


30.4 


19 L 


1, 0.6 


0.25 


186.7 


35 


3.56 


8.4 


31.5 


20 L 


1.2, 0.4 


0.5 


190.9 


49 


2.87 


4.9 


41.1 



flux characteristics is periodical in $o, what gives the 
needed assignment between measured flux and injected 
current. Thus, the coupling inductance can be deter- 
mined. In the next step, the coupling inductance was 
calculated with an algorithm described in2£, which takes 
into account geometrical and kinetic inductance as well. 
The London penetration depth A^ of the high-T c film 
acts as a fitting parameter to meet the measured cou- 
pling inductance. This A^ could be determined very 
consistently to (462 ± 11) nm for all the 20 SQUIDs. 
Now, using this known London penetration depth, the 
whole SQUID inductance could be reliably calculated in 
the same way. 

The normal resistance of the SQUIDs, Rsq was deter- 
mined from the steepness of the resistive part of VCC at 
large currents (I > 10 Isq). Now, knowing the normal 
resistance of the SQUID, its critical current Isq can be 
assessed. The direct extraction from the experimental 
VCC is difficult due to the large level of thermal noise, 
which gives strong rounding of the curves. To get access 
to the noise-free intrinsic critical current of the SQUID, 
a set of theoretical VCCs is calculated, which were ob- 
tained from an analytical solution of the SQUID equa- 
tions in the small inductance limit 28 (sec Fig. 2]). In 
comparison to this set of theoretical VCCs (current nor- 
malized to Isq vs. voltage normalized to the IsqRsq 
product) the experimental VCC is plotted. Because Rsq 
is known, the noise-free critical current Isq can be tuned 
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SQUID #9 




TABLE II: DC SQUID parameters for the bias current Imax 
given in Table I, where maximum swing AVmax is provided; 



is the flux shift of the VFC under reversal of /„ 



FIG. 5: Experimental VFC for SQUID #9. A measured flux 
shift A$ me as is determined as a distance between the first 
extremums of VFCs relative to zero flux point (I co u = 0). 



until the experimental curve fits well to any of the theo- 
retical ones. 

The maximum voltage modulation AV max and the bias 
current I max needed for that voltage swing could directly 
be taken from the measurement of the voltage-flux char- 
acteristics. 

At the bias current \ m axi which gave the maximum 
swing of the VFC (both shown in Table the next 
SQUID parameters were measured. They are presented 
in Table [TTJ The transfer functions represents the 
maximum steepness of VFC at their left (+) and right 
(-) sides. As VFC is a $o periodic function, a unique de- 
termination of the actual shift A$> max , which is given in 
the second column of Table HT1 is not, in general, possible. 
Below in Section IVB we describe in detail a reasonable 
procedure which relates A<& max with A<i> meas defined as 
the distance between the first extremums of VFCs rela- 
tive to zero flux (I co il — 0, see Fig. [5]). 

The asymmetry parameters of the SQUIDs are given 
in Table [fal 

Unlike the method of Ref. Il7|, where the asymmetry of 
the critical current was assessed indirectly from the flux 
shift of the VFC, we determined the asymmetry param- 
eters by a direct measurement of critical current and re- 
sistance of one of the two junctions in the SQUIDs. This 
got possible due to the opening of one superconducting 
arm of every SQUID loop. 

The first parameters given in Table Hill the critical cur- 
rent Ic and resistance R are average junction values for 
the SQUID as defined in Q, ©. They differ from the 
SQUID values given in Table U by a factor of 2. For a 
symmetric SQUID these are just the values of each of 
the two junctions. The quantities I 2 and R 2 are critical 
current and resistance of the second junction which re- 
mained for the measurement after the SQUID loop had 



lepresents the maximum steepness of the VFC at their left (+) 
,nd right (-) sides; AVmax / 1 R is the maximum voltage swing 
elated to the average junction resistance (see Table III). 



SQUID# group A<IW 


v£ 




AVmax 1 R 




$o 




nV/$a 


HV/Q. 


1 S 


0.357 


115.7 


-99.7 


4.531 


2 S 


0.532 


170.1 


-147.3 


4.163 


3 S 


0.675 


71.3 


-60.01 


5.234 


4 S 


0.23 


129.5 


-116 


4.855 


5 S 


0.412 


120.13 


-167.3 


4.548 


6 S 


1.243 


102.89 


-121.4 


5.236 


7 M 


1.36 


74.64 


-59.97 


2.183 


8 M 


0.079 


50.74 


-48.11 


2.123 


9 M 


1.656 


57.49 


-41.19 


1.798 


10 M 


-1.124 


52.44 


-51.26 


2.056 


11 M 


0.81 


36.57 


-37.72 


2.341 


12 M 


0.842 


45.53 


-45.08 


2.437 


13 M 


1.344 


49.85 


-37.34 


2.215 


14 M 


1.985 


52.38 


-41.92 


2.187 


15 M 


1.792 


41.95 


-34.51 


2.482 


16 L 


0.947 


22.6 


-22.29 


1.060 


17 L 


0.948 


22.58 


-25.24 


0.851 


18 L 


-0.415 


17.74 


-20.19 


0.979 


19 L 


1.455 


23.96 


-23.69 


1.181 


20 L 


2.934 


13.49 


-15.79 


0.850 



been cut. Using all these data, we obtain the bulk asym- 
metries 7 and p using equations ([1]) and (J2J), and, finally, 
with the aid of equations ([3]), the intrinsic asymmetries 
ctj and a p . So, for the first time the asymmetry parame- 
ters of the current density and resistivity of the junctions 
in high Tc DC SQUIDs were determined independently 
and directly. 

Before we discuss the influence of junction asymmetry 
on the output SQUID characteristics we should like to 
show to what extent our high Tc Josephson junctions 



obey the scaling law IcR ~ j 



1/221 



In this case the pa- 



rameters aj and a p are interrelated to each other by Eqs. 
Oi. A theoretical dependance of ay on a p (second of Eqs. 
[5]) for this scaling law together with experimental points 
is shown in Fig. [SJ As is seen from this figure most of 
the experimental points are grouped over the theoreti- 
cal curve. Therefore, we may conclude that our high Tc 
Josephson junctions are described by the aforementioned 
scaling law with a good accuracy 



IV. DISCUSSION 

A. The depth of the voltage modulation 

In general, the increase of the SQUID inductance leads 
to a decrease of the voltage swing of VFCA. This is illus- 
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TABLE III: Current and resistance asymmetry of DC SQUID. 



SQUID# 




R 


IC2 


R.2 


7 


P 




OLp 


group 


/iA 


o 
M 


. < A 
/iA 


o 
M 










1 S 


18 


7.9 


11 


9.45 


0.39 


0.16 


-0.39 


-0.16 


2 S 


11 


11.58 


g 


13.56 


0.18 


0.15 


-0.18 


-0.15 


o o 


1 7 ^ 

-L 1 .O 




7 8 


1 9 OQ 


n w 

\J.OO 


U.O -L 


n 9fiQ 


n 91 fi 

U.Z L-U 


4 S 


1 4 7^ 

J. 4:. 1 O 


7 QQ 


12.5 


o.uy 


n i ^ 

U-4-O 


u.uo 


U. -LO 


u.uo 


5 S 


12.75 


9.63 


8.5 


13.05 


0.33 


0.26 


-0.087 


-0.011 


6 S 


21 


6.8 


12 


12.77 


0.43 


0.47 


0.089 


0.039 


7 M 


15.5 


9.85 


4.2 


23.84 


0.73 


0.59 


-0.73 


-0.59 


8 M 


16.5 


7.73 


15 


8.32 


0.09 


0.07 


-0.09 


-0.07 


o i\/r 
y ivi 


on 


o i o 
y. iz 


q 
o 


1 fi 7*7 
lo. 1 1 


U.oo 


U.ol 


n 70Q 

-U. / zo 


n oi « 


10 M 


1 ^ 9^ 

-LO.ZO 


8.22 




9fS 




u.uu 


U.OL. 


u.ou 


11 M 


28.25 


5.21 


9D ^ 

ZU.O 


u.uo 


n 98 


0.14 


n 9s 

u.zo 


-0.14 


12 M 


97 ^ 


(S 9 
u.z 


24 


7 1 "3 


n i ^ 

U. -LO 


n i ^ 

U. -LO 


0.124 


0.124 


13 M 


17.5 


6.54 


5.9 


13.71 


0.76 


0.52 


-0.63 


-0.31 


14 M 


21 


7.09 


3.6 


25.74 


0.83 


0.72 


-0.564 


-0.344 


15 M 


26.5 


5.52 


9.8 


12.43 


0.63 


0.56 


-0.19 


-0.083 


16 L 


14.5 


8.68 


6.5 


13.21 


0.55 


0.34 


-0.55 


-0.34 


17 L 


12.75 


9.35 


5.3 


13.62 


0.58 


0.31 


-0.309 


0.022 


18 L 


23 


6.82 


28.5 


5.58 


-0.24 


-0.22 


0.24 


0.22 


19 L 


17.5 


7.11 


8.8 


11.46 


0.5 


0.38 


-0.286 


-0.144 


20 L 


24.5 


5.74 


6.4 


16.69 


0.74 


0.66 


-0.381 


-0.239 
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FIG. 6: The dependance of a.j on a p . The theoretical depen- 
dance was calculated from the second of Eqs. [5] 
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FIG. 7: The voltage modulation as a function of the critical 
current. The lines for L=70 pH (O), 130 pH (y), and 190 
pH (A) are calculated from Eq. [8] The experimental points 
are shown by x for S-SQUIDs, + for M-SQUIDs, and * for 
L-SQUIDs. 
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FIG. 8: The depth of the voltage modulation as a function of 
the critical current asymmetry 7. 



trated in Fig. [7] where the voltage modulation AV max /R 
is shown as a function of the junction critical current. 
As is seen from the figure the SQUIDs are grouped by 
their inductance near corresponding lines obtained from 
the expression of Enpuku22. 



AV 



I C R tt(1 + /3) 



exp — 3.57T 



2 ksTL 



(8) 



We also investigated the dependance of the maximum 
voltage modulation on the asymmetry parameters of the 
DC SQUID (see Figs. [Hand©. We have found that for 



M- and L-SQUIDs the influence of junction asymmetry 
on the voltage modulation is weaker than that for small 
inductance SQUIDs. For S SQUIDs there is appreciable 
scattering of AV ma x versus 7 and p. 



B. Flux shift 

As is known the asymmetry of the SQUID junctions 
results in a flux shift of the VFC under bias reversal. 
Here we determine the flux shift as the sum of the shifts 
of two VFCs, one for / = +I m ax and the other for / = 



7 




FIG. 9: The depth of the voltage modulation as a function of 
the resistance asymmetry p. 
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FIG. 10: The flux shift as a function of the current asymmetry 
7. Geometrically symmetric SQUIDs (a g = 0) are shown by 
black boxes, SQUIDs with a g 7^ are shown by open circles. 
Theoretical points calculated from (0 are shown by stars. 



—I maxi relative to the zero flux point (I co u — in Fig. 
[5]). Here the subscript "max" means the bias current 
which provides maximum swing of the VFC. Since the 
VFC is $0 periodic, the flux shift can be determined 
with an accuracy 2n$o, where n is an integer. Thus the 
minimum uncertainty of the determination of the flux 
shift for asymmetrical SQUIDs is 2$o- 

In an ideal case when there is no parasitic trapped 
flux in the loop, both VFCs are shifted symmetrically by 
the same amount A^> sym in opposite directions with re- 
spect to zero flux point. In this case the overall flux shift 
is calculated as a sum of both shifts, A<!> = 2A$ SJ , m . 
However, in our cases the shift of two VFCs is not sym- 



FIG. 11: The flux shift as a function of the resistance asymme- 
try p. Geometrically symmetric SQUIDs (a 9 = 0) are shown 
by black boxes, SQUIDs with a g 7^ are shown by open 
circles. Theoretical points calculated from are shown by 
stars. 

metric with respect to zero point (see Fig. [SJ. The rea- 
son for this is the parasitic flux <& p trapped in a loop. 
This flux does not change its sign under bias reversal, 
hence the flux shifts are as follows: A<I>_ = A$ sym 
A$ + = - A$ sr ,„+$ p , where A$_ is the flux shift for the 
VFC for reversed bias current (/ = — I max ), and A$ + is 
the same quantity for the unreversed VFC (/ = +I m ax)- 
Thus, the parasitic flux disappears from the overall flux 
shift A$ = A<I>_ - A$ + = 2A§ sym . The application of 
this simple formula requires precise knowledge of the cor- 
responding extremums of the VFCs shifted by a parasitic 
flux. In other words, we must know these two extremum 
points on both VFCs whose distance provides us with the 
overall flux shift. However, since VFC is $0 periodic, it 
is not possible to definitely determine these two points. 
Therefore, a single valued determination of the flux shift 
requires some additional information. 

The procedure we used for the unique determination 
of the actual flux shift, A$ act was performed in three 
steps. First, we measured the flux shift A$ meas between 
the nearest extremum points of VFCs with respect to the 
zero flux point. Second, by using expression ([7]) we cal- 
culated for a given SQUID a theoretical value of the flux 
shift A<& t /i eor . Finally, the actual flux shift A$ act was de- 
termined by subtracting from (or adding to) A$ meas the 
integer number of flux quanta n<i>o [n = 0, ±1, ±2, etc.) 
which gave the flux shift nearest to the theoretical value 

A&theor- 

The illustration of our method is presented in Table 
IIVI The values of actual flux shift from the last column 
of this table are given as A& max in the second column of 
Table El 

We want to stress here that in our method the theoret- 
ical estimations of the flux shift were used only as a guide 
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TABLE IV: Determination of actual flux shift. A$ meas is the 
flux shift measured between nearest extremums of two VFCs, 
A&theor is the flux shift calculated from Eq. [7] and A"J> ac t is 
the actual flux shift corrected for parasitic flux. 



SQUID# group 






A$ act 




$0 


$o, Eq.(7) 


$o 


1 S 


0.357 


0.374 


0.357 


2 S 


1.532 


0.203 


1.532-1=0.532 


3 S 


1.675 


0.913 


1.675-1=0.675 


4 S 


1.23 


0.157 


1.23-1=0.23 


5 S 


0.412 


0.347 


0.412 


6 S 


1.243 


0.784 


1.243 


7 M 


1.36 


1.328 


1.36 


8 M 


0.079 


0.177 


0.079 


9 M 


0.656 


1.302 


0.656+1=1.656 


10 M 


-1.124 


-1.07 


-1.124 


11 M 


0.81 


0.663 


0.81 


12 M 


0.842 


0.451 


0.842 


13 M 


1.344 


1.207 


1.344 


14 M 


0.985 


2.116 


0.985+1=1.985 


15 M 


0.792 


2.002 


0.792+1=1.792 


16 L 


0.947 


0.984 


0.947 


17 L 


0.948 


0.75 


0.948 


18 L 


-0.415 


-0.86 


-0.415 


19 L 


0.455 


1.315 


0.455+1=1.455 


20 L 


0.934 


2.8 


0.934+2=2.934 



for the determination of the actual flux shift caused by 
the predetermined asymmetry. This differs from many 
other papers on the subject (see, for example, Ref. fl7l ) 
where theoretical expression (j6]) itself has been used for 
the determination of asymmetry parameters. In this lat- 
ter case there are not any means which allow a unique 
discrimination between the unshifted curve and the curve 
shifted by a parasitic flux. 

The dependance of experimental flux shift (second col- 
umn in Table [TTJ) on current asymmetry 7 and resistance 
asymmetry p is given in Figs. llOl and fTTl Also we show on 
these graphs the theoretical points (star symbols) which 
were calculated from Eq. (J7JI by using the measured 
values of 7 and p. As is seen from these graphs, the 
equation (0 is a good approximation for the measured 
flux shift. The experimental points for the flux shift can 
be approximated by least mean square fits which are 
as follows: A$ max /<f> = 1.897 - 0.03 for Fig. [TO] and 



A^maJ^o = 2.24p + 0.12 for Fig. \TT\ 



Transfer function 



For every SQUID the transfer functions V§ were mea- 
sured (see Table [TTJ) and their dependance on different 
SQUID parameters was investigated. The dependance of 
Vgf on the loop inductance L is shown in Fig. [T3J We also 
show on this plot the transfer function for a symmetric 
shape of the VFC by using the relation V$ = 7rAU/<f>o 
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FIG. 12: Dependance of the transfer functions on SQUID 
inductance. 



with AU being taken from Table |TJ It is seen that, in 
general, the increase of inductance results in a decrease 
of both the transfer functions and the scattering for the 
transfer function values. The largest scattering is ob- 
served for S SQUIDs. Another conclusion is that as the 
SQUID inductance is increased the shape of VFC be- 
comes more symmetric: the transfer function for sym- 
metric shape becomes closer to the measured transfer 
functions V$ . As a measure of the shape asymmetry 
of VFC we introduce the quantity | — V§ | which is 
equal to zero for a symmetrical shape. The dependance 
of this quantity on SQUID inductance is shown in Fig. 
[TBI From this figure it is clearly seen that the shape of 
VFC becomes more symmetric as the SQUID inductance 
becomes larger. 

The dependance of the measured transfer functions V$ 
on (3 is shown in Fig. 1141 Here one can see a nearly ex- 
ponential decay of the transfer function with the increase 
of /?. 

The dependance of the normalized transfer functions 
V$ Qq/ RIc on current asymmetry 7 and resistance 
asymmetry p is shown in Figs. ll5l and [To1 As is seen from 
these figures, the L SQUIDs show a weak dependance on 
the parameters of asymmetry. It conforms with Fig. 1131 
where the VFCs for these SQUIDs have minimum dis- 
tortion. However, for SQUIDs with small inductance (S 
SQUIDs) the dependance on asymmetry parameters is 
much stronger. It is worth noting that the 7 and p de- 
pendances of U $ (see Figs. [15l Unj) have indications of 
a peak for S SQUIDs near the value of 0.2. A similar- 
behavior has been found for SQUIDs with L m 20 pH in 
Ref. [13 (see Fig. 6 therein). 

We also compare the dependance of measured transfer 
functions V$ on current asymmetry | aj | for geometri- 
cally symmetric SQUIDs (a a = 0) with the correspond- 
ing expressions from Ref. [l7| obtained from computer 
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FIG. 13: The deviation of VFC shape from symmetrical form. 
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FIG. 14: Dependance of the transfer functions on j3. Black 
symbols are for V^J", white symbols are for . 



simulations. The results are shown in Fig. 1171 Again, 
the strongest deviation of experimental points from the- 
oretical ones is observed for S SQUIDs (only SQUID =#=7 
is an exception). 



V. CONCLUSION 

The main purpose of the study was the experimen- 
tal investigation of the influence of the junction asym- 
metries on the SQUID output characteristics (depth of 
modulation, flux shift, transfer function). To this end we 
directly measured the current and resistance asymme- 
try of the junctions in every of 20 investigated SQUIDs. 
It turned out that the values of current and resistance 
asymmetry (a,j and a p ) were randomly distributed over 



FIG. 15: Dependance of the transfer functions on current 
asymmetry 7. Black symbols are for V£ , white symbols are 
forVT. 
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FIG. 16: Dependance of the transfer functions on resistance 
asymmetry p. Black symbols are for V£ , white symbols are 



for V a 



a chip without noticeable correlation either with SQUID 
inductance or critical current. Even for geometrically 
symmetric SQUIDs (we had 9 such SQUIDs) there was a 
significant asymmetry in critical current and resistance. 
Nevertheless we achieved definite conclusions about the 
influence of the junction asymmetry on the SQUID out- 
put characteristics. 

The dependance of the depth of modulation AV and 
the transfer functions on the junction asymmetry is 
appreciable only for low inductance SQUIDs (see Figs. 
El O HS1 [HI). For large inductances this dependance is 
rather weak which is in accordance with a small distor- 
tion of the shape of VFC for these SQUIDs (see Fig. 
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# L(pH) 




1 71.4 


18 


2 714 


11 


4 65.4 


14.75 


7 131.6 


15.5 


8 131.6 


16.5 


10 125.7 


15.25 


11 125.7 


28.25 


16 191.5 


14.5 


18 185.8 


23 




very sensitive to the junctions asymmetry. However, for 
SQUIDs with a relatively small inductance L < 70 pH 
the dependance on asymmetry is more significant, the 
shape of VFC is more distorted and the transfer func- 
tions are well above the corresponding values for large 
inductance SQUIDs. 

As was expected, the flux shift of the VFC is more 
sensitive to the junction asymmetry than the depth of 
modulation or transfer functions. The dependance of the 
flux shift on 7 and p is approximately linear and it is well 
described by the analytical model (Figs. fTTJl PTTj) . 

The results obtained in the paper are important for 
the implementation in the sensitive instruments based 
on high Tc SQUID arrays and gratings. 



I a, | 
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